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We present a formulation for the internal motion of equilibrium configurations with a rotational 
Killing vector in general relativity. As an approximation, this formulation is applicable to inves- 
tigation of the internal motion of quasi-equilibrium configurations such as binary neutron stars. 
Based on this simple formulation, a condition for the general relativistic counter rotation has been 
obtained, though in the recent work by Bonazzola, Gourgoulhon and Marck, their condition for the 
counter rotation is not enough to specify the internal velocity field. Under the condition given in 
this paper, the internal velocity field can be determined completely. Indeed, in the counter-rotating 
case, we have also derived Poisson equations for the internal velocity, which take tractable forms in 
, numerical implementation. 
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Kilometer-size interferometric gravitational wave detectors, such as LIGO §J, VIRGO § and TAMA § are now 
under construction, which are increasing our expectation of direct detection of gravitational waves. Coalescing binary 
neutron stars are among the most promising sources for such detectors. The main reasons are that (1) we can expect 
to detect the signal of coalescence of binary neutron stars about several times per year [Q, and (2) the waveform from 
coalescing binaries can be predicted with a high accuracy compared with other sources Jp],|],^| ■ Informations carried 
by gravitational waves tell us not only various physical parameters of neutron stars J^gf^ but also the cosmological 
parameters [To[] if and only if we can make a detailed comparison between the observed signal with theoretical 
prediction during the epoch of the so-called inspiraling phase where the orbital separation is much larger than the 
radius of component stars JTT[ . 

To obtain such a large amount of informations of the gravitational waves, we must understand the theoretical 
O ■ mechanism of merging and construct templates for signals a priori. When the orbital separation of binary neutron 
stars is < 10GM/c 2 , where M is the total mass of binary neutron stars and c is the light velocity, they move 
approximately in circular orbits because the timescale of the energy loss due to gravitational radiation tew is much 
longer than the orbital period P as 



I. INTRODUCTION 



taw dc 2 

— *Htqgm) l^J' (L1) 

where fj, and d are the reduced mass and the separation of binary neutron stars. (As for higher order corrections, see |lj] 
for instance.) Thus, binary neutron stars evolve adiabatically radiating gravitational waves. As the orbital separation 
approaches QGM/c 2 , however, they cannot be described by the circular orbital motion because of instabilities due 
to the general relativistic gravity and/or the tidal field. As a consequence of such instabilities, binary neutron stars 
will plunge to merge. Thus, it is expected that the property of the signal of gravitational waves changes drastically 
around this transition epoch. This means that gravitational waves emitted at this transition region may bring us an 
important information about the internal structure of neutron stars, since the location where the instability occurs, 
the so-called innermost stable circular orbit (ISCO) significantly depends on the equation of state of neutron stars 

A strategy to search the ISCO accurately is as follows: Close binary neutron stars have two dominant time scale s, th e 
orbital period and the time scale of the energy loss. The former is much shorter than the latter according to Eq.( [l.l| ). 
From this physical point of view, we may consider that binary neutron stars evolve in the quasi-stationary manner, 
so that we can take the following procedure; first, we neglect the effect due to gravitational radiation and construct 
equilibrium configuration, and then the radiation reaction is taken into account as a correction to the equilibrium 
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configuration. However, how to separate the stationary part from the non-stationary part is one of the most important 
but difficult problems in general relativity. For instance, Detweiler proposed that a stationary solution of the Einstein 
equation with standing gravitational waves will be constructed by adding incoming waves from infinity, and may be 
a valuable approximation to physically realistic solutions Jlq| . However, these solutions are not asymptotically flat 
because the total energy of gravitational waves inside a radius r grows linearly with r. Thus, the boundary condition 
for them is significantly different from that of physical solutions. 

In the numerical approach, Wilson and Mathews |17| proposed a semi-relativistic approximate method in order to 
obtain the equilibrium configurations of binary neutron stars just before merging. In their method, they assumed the 
spatially conformal flat metric for binary neutron stars. Practically, following their approach, numerical calculations 
have been done by some people [fl8|-p0| . However binary neutron stars are genuinely anisotropic, so that the assumption 
of the spatially conformal flatness should break down at some level. As a result, a lot of debates are continuing [pl|-p4|. 

In the post-Newtonian approximation, the metric and the material quantities are expanded with respect to c _1 
assuming the slow motion and weak gravity p5j-|29[] . In this approximation, we can identify the radiation reaction 
terms which begin at the 2.5PN order |2S(] . Thus, it is possible to construct the equilibrium configuration of binary 
neutron stars in the 2PN approximation |3p| , |3l|l . Asada and Shibata presented a formalism to obtain equilibrium 
configurations of uniformly rotating fluid at the 2PN order pl[ . In this formalism, the equations to be solved are 
reduced to Poisson equations. This formulation is applicable to construction of non-axisymmetric uniformly rotating 
equilibrium configurations which include synchronized binary neutron stars and the Jacobi ellipsoid. 

Nevertheless, owing to Kochanek's investigation, the viscosity inside close binary neutron stars may be significantly 
small, so that synchronization may not be a good approximation |32[ . Rather, irrotation (counter rotation) may be 
preferred for close binary neutron stars. Recently, Bonazzola, Gourgoulhon and Marck presented a formulation of the 
general relativistic Euler equation for perfect fluid hydrodynamics under the maximal slice ]33| ] . Especially, in the case 
of the counter rotation, they gave Poisson equations to determine internal velocity fields, where they put only one 
condition on some scalar so that the counter rotation could be chosen from various kinds of internal motions allowed 
by the equation of motion. However, in order to specify the internal spatial velocity field which can be expressed 
as an essentially three vector, we must adopt three conditions. Indeed, although one of their equations (Eq.(47) in 
p3|) is the Poisson equation for the internal velocity, it turns out to be the equation projected perpendicularly to 
the internal velocity. This will be described in detail in section 7. Therefore, only their equation is insufficient to 
determine the velocity field. The main purpose of this paper is to investigate carefully the equations to determine 
the internal velocity field with respect to the orbital motion, and present the equations to determine completely the 
internal velocity field in the case of the general relativistic counter rotation. In section 2, we assume the rotational 
Killing vector and introduce basic quantities. We consider the conservation law with the above point in mind in section 
3. In section 4, the counter-rotating case is considered. In section 5, some basic equations are derived, and based 
on this formulation, the condition for the counter rotation is presented. Next, equations to determine completely the 
internal velocity field are derived in section 6. In Section 7, we discuss the relation with this formulation with the 
previous work. Section 8 summarizes the conclusions. Appendix shows the relation between the counter rotation and 
the irrotation. 

We use the units of c = G = 1 in this paper. Greek and Latin indices take 0, 1, 2, 3 and 1,2,3, respectively. 

II. KILLING VECTOR AND TWO FRAMES 
A. Killing vector 

Motivated by the idea discussed in the Introduction, let us consider the equilibrium configurations by assuming the 
Killing vector such as 

l a =t a + nm a , (2.1) 

dt 



where t a = {4i) a and m a = {-§r) a ■ If one uses the Cartesian coordinate, m a may be expressed as 



The Killing equation is 



where C denotes the Lie derivative. 



Cig a /3 = VJ P + Vpl a = 0, (2.3) 



We consider a foliation of the spacetime described by the hypersurface (E t ) whose unit normal vector is 

n^ = (-N,0), (2.4) 
1 N* 

n» = (jJ,- W ), (2-5) 

where N and N l are the lapse function and the shift vector respectively J|4||35|]. The line element is expressed as 

ds 2 = -~N 2 dt 2 + g ij (dx i + N l dt)(dxi + N j dt), (2.6) 

where our sign convention for the shift vector is opposite to that of Bonazzola et.al |33|. We consider the internal 
motion which may deviate from the orbital motion (co-rotation). Therefore, it is convenient to introduce two frames 
for the purpose of the detailed investigation |33| ; one is the non-rotating frame and the other is the co-rotating frame. 

(1) non-rotating frame 

At first, we consider the non-rotating frame as follows: The temporal vector t a is related with the unit normal 
vector n a as 

t a =Nn a + N a . (2.7) 
Then, in this frame, we rewrite the Killing vector l a as 

r = Nn a + B a , (2.8) 

where we defined B a as 

B a =N a + VLm a . (2.9) 

We define the extrinsic curvature as 

K af 3 = -(P«)S(^.)^V M n„, (2.10) 
where the projection tensor (P n )a was defined as 

(P n )Z = S£ + n a nf>. (2.11) 

Then, we obtain |3(| 

V a »/s = -K af3 - n a D/3 In N, (2.12) 

where D a denotes the spatial covariant derivative defined on E t and we used n^dx* 1 = —Ndt. 

(2) co-rotating frame 

The decomposition with respect to l a is useful to describe the internal motion relative to the co-rotation along l a . 
Let us introduce the unit tangent vector along l a as 

„« = e -*ja (2.13) 

and the projection tensor (P v )a as 

(Pv)Z = SZ + v a v0. (2.14) 

Thus, in the co-r otating frame, we use the projection tensor (P v )a> m place of (P n )a m the non-rotating frame. 
From Eq.( 2.13| ), v a v a = —1 becomes 

e 2 * = N 2 - B a B a . (2.15) 

From Eq.Q2.13D, wc obtain @ 

Vpv a = uj a/ 3 - w^V Q $, (2-16) 
where using the anti-symmctrization [ ] wc defined 

^ a p = -(PvY a (PvY^ [ ,v v] . (2.17) 
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Here, introducing 



we define the vorticity of v" as 



e°0T = v^ al3 \ (2.18) 



uj a = --e a ^LU^, (2.19) 



where we used 

w M " 7 = (PvUPv)p ~ (PvZ(p v r p . (2.20) 

III. CONSERVATION LAW 
A. adiabatic flow 

Since we consider the perfect fluid in this paper, the stress energy tensor is written as 

T" v = {e + P)u^u u +Pg»". (3.1) 
We decompose the four velocity in terms of v a as 

«/* = r(w" + V"), (3.2) 

where we defined T and V^ L as 

r = (3.3) 
v^* = r-^P^^u". (3.4) 

Then, from u^u^ = — 1, we obtain 

r = (i - V^V 1 )- 1 ' 2 . (3.5) 

Next, for the adiabatic flow, the first law of the Thermodynamics leads to 

V Q e V Q n 
e + P n 

where we assumed the barotropic equation of state, e = e(n) and P = P(n). Instead of the specific enthalpy 



(3.6) 



e + P 

H=^—, (3.7) 
rriBn 

where mg is the baryon rest mass, we use the logarithmic enthalpy as 

h = bxH. (3.8) 

Eq.(|3.(j) is rewritten as 

V^=^- (3.9) 



4 



B. conservation law 



Here, we consider the conservation law J33|. From its decomposition with respect to v u , we obtain 

V M ^+V^V M (ln(r 2 ri) + $ + /i), (3.10) 

and 

V M (FV a ) + 2oj atl V fi + V Q $ + T- 2 V a h 

+V a V» V M (in (Tn) - = 0. (3.11) 



Moreover, by decomposing Eq.(3.11) with respect to V a , we obtain 

7%(lnr + $ + /i) = 0, 



(3.12) 



which can be taken as the expression for a general relativistic version of the Bernoulli theorem. Namely, lnT + <f> + h 
is constant along the internal flow V a . Then, Eq.(3.10) becomes 



V V" + ^V M ln(nr) = 0, 



(3.13) 



which is exactly the equation for the baryon number conservation. Furthermore, we obtain the general relativistic 
analogue of the Euler's equation as 



(3.14) 



where we used Eq.(3.12). 

In order to obtain simple expressions later, it will be useful to introduce the vorticity of the internal motion as 



Using the internal vorticity, we obtain the relation 

v M u Q = r~ 2 v a for + 2e afiV n l ,v v , 



where we used Eq.(3.5). Using this relation, Eq.(3.14) becomes 

V Q (lnr + $ + h) + T 2 (2l a ^(n^ + uJV v + V^fPy^V 1 ^) = o, 
where we defined the projection operator with respect to the internal velocity as 



{Pv)i = € 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



IV. COUNTER-ROTATION 

In order to simplify the expression, we introduce the scalar Q as 

g = inr + $ + /i. 

This can be re-expressed as 

Then, Eqs.( 3.12] ) and ( 3.17j ) become respectively 

v^v^g = o, 



(4.1) 



(4.2) 



(4.3) 



and 



2e a ^(n" + wjv + F^F^Py^v^ + r~ 2 v Q g = 0. 



It should be noted that we do not need to solve Eq.(4.3), since it is obtained by contracting Eq.(4.4) with V a . 
In the Newtonian limit, the counter rotation means 

fin + — > o. 

Therefore, a simple condition for the counter rotation would be 

ft M + = 0. 



In this case, however, Eq.(4.4) becomes 



r 2 v^(p v )»v„<s> + v a g = o, 



(4.4) 

(4.5) 
(4.6) 
(4.7) 



for which such a Q does not always exist. 

It seems difficult to find o ut c onditions for the internal motion in the counter-rotating case from Eq.(4.4). Therefore, 
it is desirable to change Eq.(4.4) into the form easier to investigate. Indeed, we can reformulate the equation as follows: 
First, we should note the relation 



where we defined the renormalized vorticity as 

tt a = 

Then, we obtain 

I, 



-e*e"^V M (e-*K) 



Using this relation, we rewrite Eq.(4.4) in the simple form as 

ea^n" + ^ l )v» + ^r- 2 v a g = o. 



Owing to the simple form of Eq. ( 4.11 ) , we can define the counter rotation by requiring Q." such as 

n" + w" = 0. 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4-12) 



In t he N ewtonian limit, this coincides with the condition for the Newtonian counter rotation. Under Eq.( 4.1 j ), 
Eq.( 4.1l| ) means simply 



Q = constant, 



(4.13) 



which satisfies a general relativistic version of the Bernoulli theorem (3.12). It is noteworthy that the condition ( 4. 12) ) 
satisfies the irrotation condition. (For instance, see JItJ for the definition of the vorticity). This is shown in the 
appendix. 



V. BASIC EQUATIONS FOR INTERNAL VELOCITY FIELD 
A. decomposition of V a 



In the previous section, we h ave presented Eq.( 4.12j ), the condition for the general relativistic version of the counter 
rotation, by the use of Eq.( [4.4| ) which is an algebraic equation for the vorticity fl a . For the numerical implementation, 
it is desirable to derive the equations for the internal velocity field in the case of the counter rotation. First, with 
respect to S t , we decompose the internal velocity V a into the following form 



G 



where we defined 



V( = V\\n a , 
Vf| = -npVl 3 . 



(5.1) 

(5.2) 
(5.3) 
(5.4) 



Since the spatial vector V a expresses only three degrees of freedom in the four velocity, Vj| is related with V". 
Indeed, from v a V a = 0, we obtain 



For the sake of the further decomposition, we introduce the spatial alternating tensor as 



This new alternating tensor is related with t a ^ as 

In turn, with respect to n a , we decompose the covariant derivative of the internal velocity as 



V a Vp = D a V^f3 + -±[B, V ± ] p + (naKfr - n K ai )V2 - n a n p VlD n InN 
Vj, - V\\{K a0 +n a D \nN), 

where we defined the Lie bracket as 

[B, VA_] a = B t5 VpVZ - V[ VpB a . 
Here, we used the following relation for the vector X a on S t 



/V,jl a = n a X ts D f) \D.N - XpK al1 - — [B,X} a , 



a/3 



1 



which is derived from 



C,X a = 0. 



Therefore, by the frequent use of Eq.fljUCl) in Eq.Q, we obtain 
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ND^Vx v + 5 M V1U„ - 2if„ TT - V\\D V InN + D V V\\ } 



NVxr-VwBp) - —V_ Lll B v B' T D a $ 



-e-^n a e x ^B x [D lx V Xu + V ,/) \> 



where we defined the Lie bracket on £ t as 

[B,V±]± u = Bi 1 D lt V± v -V£D lt B u . 
From Eqs. Q2.17D , d 2 - 19 D and 10), we obtain 

c^" = ^e-^e^^ND^ + 2B^D V N + 2K°B a B^j + i e " 2 *n Q B„?» v D „.B„. 
From Eqs. Q5.12p and ( |5. 14j ) , Eq.( |4.12| ) becomes 

e a ^ [nD„Vx„ + B^~[B, V±]± v - 2K va VZ - V\\D V InN + D V V^ 
+ (nV ±ii - V\\B^D V ^> - ^Vx^B'D^ 
+ e~* (NDpBv + 2B^D U N + 2K°B a B u ) 
+n a e x ^B x (D^Vx „ + Vl + e-*D^) = 0. 



(5.5) 
(5.6) 
(5.7) 



(5.8) 
(5.9) 

(5.10) 
(5.11) 



(5.12) 
(5.13) 
(5.14) 



(5.15) 



This equation has been already decomposed manifestly with respect to n a . Therefore, using the spatial indices, 
equations to be solved can be rewritten as 



1 



NDjV± fc + Bj (— — [B, V±]±k - 2K M V[ - Vu D k In N + D k Vu 



N 



NV ±j - V\\B D )D k § - -V ±J B k B l D^ 



-e~* (NDjBk + 2BjD k N + 2K l 3 B t B k 



= 0, 



and 



tijk 



B^DjVxk + V ±j D k <f> + e-*Dj-B fc ) = 0. 



(5.16) 



(5.17) 



It should be noted that we need not solve the latter equation, since it can be derived from the projection of the former 
onto B. L . 

we obtain 



By taking the contraction of Eq.(5. 



y a V a = D a Vl + V?D a hxN - ^-D a V\\ - V\\K, 

where we used the following identity for arbitrary vectors E a and F a on S t 

n a [E,F] a = 0. 



Thus, with the spatial indices, Eq.(3.13) becomes 



Mi 



DiVi + VIDi In (NnT) - -j^B l D l In (V^nT) -V\\K = Q. 



(5.18) 



(5.19) 



(5.20) 



Since we have derived the basic equations ( 5.16| ) and ( 5.20Q to be solved, all we must do is to investigate the 
boundary condition for these equations. Here, we concentrate our attention on conditions on the surface of the 
matter. By multiplying Eq.(5.20) with n and evaluating it on the surface of the matter, we obtain 



B l B^Vj_ : 
N 2 ^ 



(D in ) 



surface 



(5.21) 



where we used the fact that n must vanish on the surface. This is the boundary condition for Eqs.(5.16) and (5.20|). 



VI. POISSON EQUATIONS FOR THE INTERNAL VELOCITY FIELD 



In order to reduce the problem of solving Eqs.(5.16) and ( |5.2C| ) to that of solving the Poisson equations, we 
decompose V]_ as follows; 



V[ = ? ]k D 3 A k + Dty, 
where in order to fix gauge freedoms in this decomposition, we adopt the Coulomb gauge 

D,A l = 0. 

Applying this decomposition to Eq.( 5.16| ), we obtain 
I), !)■' A' = 3 R\A k 



ijk 



fc V|| 



Ml 



1 



V Xj ^B 3 )D k <5> - —V X3 B k B l D^ 



N 2 



DjB k + jjBjDkN + —KjBiBh 



N 



(6.1) 
(6.2) 



(6.3) 
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where 3 R z k denotes the Ricci tensor defined on E t and we used 



Similarly, from Eq.( |5.2C| ), we obtain 



e ijk DjV X k = -DjD 1 A 1 + ' 6 R l k A 



3ni \k 



DjD^tfj = -VlDi In (NnT) + j^B l D t In (V l{ nT) + V\\K. 



(6.4) 



(6.5) 



VII. DISCUSSION 

Bonazzola et.al. fl33f| have recently developed the formulation for the internal motion of quasi-equilibrium configu- 
rations. However, they put only one condition for a scalar, so that the counter rotation could be chosen from various 
kinds of internal motions. Namely, in our notation, their condition is expressed as 



Q = constant, 



(7.1) 



which is not enough for full specifying the internal vorticity. We can under stan d the reason by counting the degrees 
of freedom: Although the spatial velocity has three degrees of freedom, Eq.( |7.l| ) fixes only one degree of freedom. 

We can also realize the situation in the following concrete form: The equation (Eq.(50) in pq|), which was claimed 
to determine the internal velocity field, is rewritten as 



(7.2) 



which is obtained by contracting Eq.(4.4) for Q — constant with e' x " V^. Since this manipulation is the projection 
perpendicular to V a , the resultant equation (7.2) shows clearly the projection, though in the original form (Eq.(50) 
in p3|) the projection of fi Q is separated into two parts, one in the left hand side and the other in the right hand side. 
Therefore, this equation is not enough to determine fu lly th e internal velocity field. On the other hand, in this paper, 
we have defined the counter rotation by requiring Eq.(4.12), so that the full determination of the internal motion can 
be achieved by Eqs.(6.3) and (6.5). 



VIII. CONCLUSION 



In this paper, we have presented a formulation for the internal motion of equilibrium configurations with a rotational 
Killing vector in general relativity As an approximation, this formulation is applicable to investigation of the internal 
motion of quasi-stationary co nfigur ations such as binary neutron stars. In particular, by careful treatment, we have 



derived the simple equation ( 4.11 ) to determine the internal velo city relative to the orbital motion. Thus, using 



this simple form of the equation, we could find easily a condition (4.12) for the general relativistic counter rotation. 
This condition determines completely the internal velocity field. In this counter-rotating case, we have also presented 
Poisson equations for the internal velocity, which take tractable forms in numerical implementation. For instance, an 
iterative procedure has been proposed [|33fl . Thus, the formalism presented here will be useful to investigate quasi- 
equilibrium configurations for counter-rotating binary neutron stars and/or the Riemann ellipsoid. In addition, it is 
noteworthy that slice conditions are not specified in the present paper. 

We focused our attention on the conservation law, which gives us the equation of motion. In order to construct 
quasi-equilibrium configurations, we must solve the Einstein equation. The post-Newtonian approximation is useful 
to find out the stationary parts in the Einstein equation. In the counter-rotating case, the equations to determine the 
metric will be derived in the second post-Newtonian approximation, in the similar manner in the synchronized case 
PH . This will be done in the future. 
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APPENDIX A: COUNTER ROTATION AND IRROTATION 



Here, we shall show that the counter rotation (4.12) means the irrotation (vanishing vorticity), in the spacetime 
with the rotational Killing vector. The vorticity is defined as (for instance, J37J) 



For the fluid with the barotropic equation of state, the vorticity is rewritten as 

= —jj57[ a (Hup\), 

where we used the projection of the conservation law with respect to u a p3[ | 

u^V^u a + V a h + u a u^\7„h = 0. 
From Eqs. ( |2.19D and ( g]| , we obtain 

le% pa Vo(e-*(V" + v°)). 



From Eqs. (3. 8), (4.1) and (4.13), we can take e as 



e~* = CTH, 



where C is some constant. Then, Eq.(A4) becomes 



Inserting Eq.(A2) into Eq. ( |A6| ) , we obtain 



Thus, we find 

e^ a f3~ ( (& + uj^) = -2CHe*(v a nfa + v Sl ia + v 7 n afj ). 
Since the left hand side vanishes in the counter rotation, we obtain 

W ct O ( g 7 + Vpil^a + = 0. 

Contracting u^v 1 with this, we find 

n ai v 7 = o. 

Hence, contracting v 1 with Eq.(pi9|), we obtain 

n a/3 = o. 
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